|ittp: / / xxx.lanl.gov / e-print / physics / 0205035 



Damping of electromagnetic waves in low-collision 

electron-ion plasmas 

V. N. Soshnikov Q 

Plasma Physics Dept., 
All-Russian Institute of Scientific and Technical Information 
of the Russian Academy of Sciences 
(VINITI, Usievitcha 20, 125315 Moscow, Russia) 

Abstract 

Using previously developed method [jl], |2| of two-dimensional Laplace transform we 
obtain the characteristic equations k(u) for electromagnetic waves in low-collision fully 
ionized plasma of a plane geometry. We apply here a new, different from the one used 
in |l], U , iteration procedure of taking into account the Coulomb collisions. The kinemat- 
ical waves are collisionally damping in the same extent as electromagnetic waves. Despite 
the different from [|[ appearance of the dispersion (poles) equation, the obtained decre- 
ments for fast and slow wave modes coincide with results obtained in if one neglects 
the terms of higher orders in t^/c 2 , [v x and c are electron and light velocities). We 
point out how one can determine mutually dependent boundary conditions allowing to 
eliminate simultaneously both the backward and kinematical waves for transversal as 
well as for longitudinal oscillations. 
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1 Introduction 

Propagation of electromagnetic waves in low-collision fully ionized plasma is described by 
asymptotic solution of the coupled kinetic and Maxwell equations. The trivial fact is known: 
an exponential solution exp(ikx — iu)t), proposed by L. Landau in 1946 (in the simplest case 
of a plane geometry) with complex oo, is not a solution of either the Vlasov equations (for 
longitudinal plasma waves), nor the equations for transversal waves. Nevertheless in the 
available literature one usually admits that namely Landau solution is true but the above 
mentioned equations must be correspondingly corrected by additional terms according to 
Landau rules of passing around poles in calculation of logarithmically divergent integrals 
appearing at substitution of the solution into the primary equations. 

The proposed method of two-dimensional Laplace transformation combined with Vlasov 
prescription of calculating divergent integrals in the sense of principal value allows one to 
obtain very simply asymptotical solutions of the original equations. 

In this work we briefly describe both the techniques and the results of the proposed in 
new iteration procedure. Following to this method one replaces Laplace image Q plV2 {v e ) of 
the Coulomb collision term Q(v e ,x,t) by the term 
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where f P ll 2 (v e ) is Laplace image of the perturbation of electron distribution function, Q PlP2 (v e ) 
is Laplace image of the collision term, calculated in the null-iteration approximation using 
Laplace image f° of the collisionless approximation of /W 



\e\ dft\v) E P1P2 , ««/W(tO 
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In our boundary problem we can suppose = 0. 

The method of subtraction of unphysical backward field waves suggested in at non-zero 
f^f and fj£ does not yet define single- valued dependence of ff?> and ffg on v e . As we said 
before, the boundary conditions are not independent and the given boundary electrical field 
E(0,t) defines the boundary function fi(v,0,t). Such an interrelation of these quantities 
can be determined through the natural condition of absence of kinematical waves. In the case 
E(0,t) = E cos(ut), p\ = ±icu, the general expression for f^(v) is 

f (f [S) = ^L + ^L, ( 3) 

v P\ +IU) pi — 10J 

where symbol * means complex conjugation. 

By equating amplitudes of the kinematical waves to zero one obtains linear integral equa- 
tions for determination of a(v) and a*(v): 



df^ E 



+ v x f^(v)=0 (4) 

P2 = TPl/Vx 



Pi = ±icu ' '" 
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where E PlP2 (see Eq.(37) in 0) contains integrals of the type 

f$u x u z d 3 u 
Pi + u x p 2 

These equations define uniquely the dependence a(v) on v. But determined in this way /, 
can not be used to eliminate unphysical backward field waves in E(x,t), as was supposed 
in [|J. To this end one must use the boundary condition for F Pl (that is Laplace transform of 
dE(x,t)/dx\ x=0 ). 

In this way function f^> has the form 



(e), 

pi 



pi 
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where rj(v) is some complicated function of v. Assuming that factor df^ /dv z is the main in 
the dependence of f p f{v e ) on v we can use Eq.(H) with replacement r]{v) — > rj{v) for rough 
estimates. Then, in the expression for collision term, f^> in Eq.(Q) can be approximately 
omitted (both terms in Eq.([|) have the same structure in v). 

Analogous considerations in the case of longitudinal waves lead to the determination of 
but there are no other free boundary conditions to eliminate the backward waves in E(x,t). 
This fact leads to the inevitable conclusion that the normal boundary component of electrical 
field E(0,t) = Eocos(u>t) is broken at the plasma boundary due to the surface charge. This 
plasma boundary field can be found with the proportional changing of f p f(v) and E — > E' Q , 
E Pl — > E' in linear equation of the type (§) without changing field amplitude E in E P1 in 
equation for E piP2 (see Eq.(37) in @). 
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2 Collisional damping of electromagnetic waves 



Characteristic equation which is an equation for double poles pi, p 2 of Laplace images of 
electrical field E(x,t) and distribution function fi(v e ,x,t) has been obtained in 0. For 
£"(0, t) ~ E exp(u;t) it has the following form 



G(PUP2) = (pi -iuj) 



d 3 v (df^ m e Q° pip2 (v) 



c 2 J pi + v x p 2 \ dv z \e\ E. 



PlP2 



0, (7) 



where Q° PlP2 (v) ~ E PlP2 (see §, Eqs.(27)-(29)). 

The integrals are defined according to Vlasov principal value prescription. The residue of 
Eq.([7D at pole p\ = iu defines poles in p 2 . 

In the case of procedure (1) the pole equation differs in form from the pole equation in [0: 



G(pi,p 2 ) = (Pi - iw) 



p\ uj 2 lPi rdffi 



where 



P2 9 1 > 
c z c z 
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dv z pi + v x p 2 - $ PlP2 (v) 



fi(v,p 1 ,p 2 ) 

does not contain the value E PlP2 ; fi(v,p\,p 2 ) is defined by Eq.(||). 
In analogy with M, one uses approximation 
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where coefficient 1/2 appears owing to difference in normalization of distribution functions 
taken in intervals (—00, 00) and (0, 00) and 



V 0x = V u x • 
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Then one obtains 
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where 



and 
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®t 1P2 { v 0w,V y ,V z ) = % ip2 (v 0x ,v y ,v z ) ± $ PlP2 (-wox,^,w 2 ) 



The replacement $ plP2 (f x , v y , v z ) — > $ pl p a (uosj ^yj ^z); ^ — ► u cL — kT/m e is made after taking 
derivatives in v x , v y , and v z in the differential operator Q plP2 {v); integrals in dv y and dv z can 
be approximately estimated by their mean values: 



00 j —00 
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-v z dv z ~ N y N z Q(v 0x , v 0y , v 0z ) (13) 



with evident normalization constants N v and iV*. 



After elementary transformations one obtains dispersion (poles) equation in the form 
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where 
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and further as 
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This equation is an analogue of the characteristic equation (30) in 
For both electron modes 
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we obtain from Eq. 
Eqs.([lJ), (UJ): 



at |<j( 1,2 )| <C 1 and neglecting terms with higher orders in Vq x /c 2 in 
<5« = ± ^Lu 2 (21) 
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that coincides with corresponding expressions for 5*- 1 ' 2 ' in [[J, in spite of differences of charac- 
teristic equations. 

Let us emphasize here the sharp increase of the dissipative collisional absorption of electro- 



magnetic waves proportional to 1/y 1 — uj 2 l /uj 2 at uj — > ujl + with the dominating collisionless 
non-dissipative reflective evanescence of waves at u < u>l- 



3 Conclusions 

Coincidence of collisional damping decrements for the two variants of iteration process is an 
evidence of the proposed calculation method correctness. The approximation (|1|) is more 
preferable than one used in since it has more evident physical sense. The difference of 

both iteration procedures appears at large values of v^/c 2 , however in this case there is also 
growing contribution of relativistic corrections to the original equations. 

The requirements of absence of both unphysical backward (divergent at x — ► oo) field 
waves and kinematical waves smearing electron distribution function f±(v,x,t) in v x lead at 
the given boundary field E(0, t) to the determination of the boundary distribution function 
fi e \v,0,t) and the solution f{ e \v,x,t) with single- valued dependence on v, x, and t. 
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